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Abstract 

We study Dick quark-antiquark potential (up to a color factor) Vo(r) = = ^ JL + 
9f y^ ^li ln[exp(2mr) — 1] in the heavy meson sector. This potential emerges from 
an effective dilaton-gluon coupling inspired from string theory and proves to be lin- 
early rising at large distances which ensures quark confinement. The semi-relativistic 
wave equation which appears in the theory of relativistic quark-antiquark bound states 
is used. This equation is cast into a constituent second order Schrodinger-like equa- 
tion with the inclusion of relativistic corrections up to order (v/c) 2 in the quark 
speeds. The resulting equation is solved for Dick potential via the Shifted-/ expansion 
technique (SLET). The obtained results show that the spin-averaged energy levels of 
heavy-mesons are well explained and agree with other potential models or QCD sum 
rules predictions. Moreover, as a by-product, our analysis assign to the dilaton a mass 
around 56.9 MeV lying within the range of many theoretical scenario of dilatonic mass. 
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1 Introduction. 



The problem of calculating spin-averaged energy levels heavy-mesons is a very old 
subject, but still an important theme in the existing literature. At present muchexperi- 
mental material on the masses of ground and excited states of heavy quarkonia has been 
accumulated. By comparing theoretical predictions with experimental data, one can ob- 
tain information on the form of the potential of the quark-antiquark interaction. This 
information illuminates the most important features of quantum chromodynamics (QCD) 
dynamics: the asymptotic freedom and confinement. At present it is not possible to ob- 
tain the potential of the interquark interaction in the whole range of distances from the 
basic principles of QCD. Consequently, the corresponding interquark potential has to be 
determined phenomenologically. Therefore, to explain the meson masses within the exper- 
imental limits, serious attempts to build a quantitative potential model for quarks were 
made by many authors. The number of papers on quark potential model is so large that 
one cannot quote even a substantial fraction of them, but simply call attention to review 
[1], which contains many references on the subject. 

To understand confinement in gauge theories, much attention has been focused on 
quark models based on QCD, because of their success in understanding the spectroscopy 
of mesons, especially the new ones. At short distances, QCD suggests a Coulomb-type 
potential — ^r, where a s is the quark-gluon fine-structure constant. At long distances one 
expects a confining potential. Lattice gauge theory [2], and string models [3] lead one to 
expect a linear confining potential. Now it is firmly established that the combination of 
linear confining potential, plus coulombic-type short distance potential, plus one-gluon- 
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exchange forces provides a good fit to meson mass spectra [4]. However, one would be 
hard pressed to say that they are well understood. Too many of the QQ states predicted 
by the quark potential model have yet to be seen-only the ground state S, and P wave 
multiplets are filled. Therefore, until confinement is better understood, more states have 
been sited, and their properties measured, we cannot say that the subject is completed. 

Recently, it has been pointed out by Dick [5,6,7] that an interquark potential is con- 
structed through the inclusion of a scalar field (dilaton) in gauge theories. This potential 
proves to be rising at large distances which ensures quark confinement, and at the same 
time it incorporates the asymptotic freedom at small distances. Therefore, the aim of this 
paper is to dedicate more efforts to understand this new confinement generating mech- 
anism through the investigation of the phenomenological application of Dick interquark 
potential Vb(r) in the heavy mesons sector. This problem will be approached as in a pre- 
vious work [8]. Therein, it has been demonstrated that the shifted-^ expansion (SLET), 
where I is the angular momentum, provides a powerful, systematic and analytic technique 
for determining the bound states of the semi-relativistic wave equation consisting of two 
quarks of masses mi, iri2, and total energy M in any spherically symmetric potential, even 
one which has no small coupling constant parameter. It simply consists of using 1// as a 
pseudo-perturbation parameter, where 1 = 1-/3, and j3 is a suitable shift. This shift is 
vital for it removes the poles that would emerge, at lowest orbital states with I = 0, in 
our proposed expansion below. This method yields very accurate and rapidly converging 
energy eigenvalue series. It also handles highly excited states which pose roblems for vari- 
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ational methods [9]. Moreover, relativistic corrections are included in a consistent way. 
In this spirit, this paper is organized as follows: Section 2 recalls the conversion of the 
semi-relativistic equation into an equivalent Schrodinger- type equation, and there the 
Shifted-^ expansion technique (SLET) for this equation with any spherically symmetric 
potential is introduced. Section 3 is devoted to describe the phenomenological application 
of Dick quark-antiquark interaction potential. We present, in this section, our numerical 
results of the spin-averaged energy levels of charmonium, bottomonium and be families, 
in connection to the dilaton mass, and then we draw our conclusion. 

2 SLET for the semi-relativistic wave equation with any 
spherically symmetric potential 



There are different methods to calculate the bound state energies of the semi-relativistic 
wave equation [10], and references therein. However, the method we present here is chosen 
for different reasons; ease of implementation, high accuracy, and substantial computation 
time reduction is achieved. An expansion in the powers of (v/c) 2 up to two terms in 
the semi-relativistic equation which is a combination of relativistic kinematics with some 
static interaction potential yields [8]: 

Id 2 1(1 + 1) , . , E ni V(r ' 



+ \ o + 7(0 + 



Rni(r) = (^+E ne )R ni (r), (1) 



2/j, dr 2 2/ir 2 7] 

where 7(r) = V(r) — V 2 (r)/2t], E n £ = M — m\ -ra2, /i = mim2/(mi +1112) is the reduced 
mass, and rj = v//J> 2 , where v is a useful parameter defined as v = m\m\i '(m\ + m^)- In 
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this reduction formalism the spherically symmetric potential V(r) which represents the 
interaction between the two particles remains unspecified. 



With the shifted angular momentum I = 1 + (3, Eq.(l) becomes 
~l + (2/3 + l)/r+/3(/? + l)/P , 7 (r) , E ne V(r) 



1 d 2 R n£ (r) p 
2/i cir 2 



2fir 2 + Q + Qq 



(2) 

where n in this paper is the radial quantum number, and Q is a constant that scales the 
potential V(r) at large-£ limit, and is set for any specific choice of I and n, equal to P at 
the end of the calculations. 

Following our earlier work [8], we present here the relevant formulas obtained in the 
SLET framework for semi-relativistic motion of a particle bound in radially symmetric 
potential. The calculations, and the results are very simple. The leading-order binding 
energy is given as: 

E o = V(r )- V + J V 2 + ^. (3) 
V » r o 

Expanding all the quantities in powers of 1/7 as described in Ref. [8], one finally gets the 
second-and third-order corrections E2, and £3 of the energy E n £ as: 

E 2 = 9m (4) 

where and «( 2 ) appearing as a correction to the leading order of the energy expression 
are given in the appendix of Ref. [11]. 



Collecting these terms, and carrying out the mathematics immediately gives an ex- 
pression for the energy eigenvalues, that is 

a (i) 



E n e — Eq + 



+ 



a 



(2) 



+ O 



(6) 



where I = I — (3, (3 is chosen so that the next contribution to the leading term in the 
energy eigenvalue series to vanish, i.e., E\ = 0, which implies that (see Ref.8) 



with 



and Q satisfies 



ijj = — 



/3 = -l/2-/i(n + l/2)w, 



3 + r V"(r )/V'(r Q ) - fir^V (r Q y / (Qrj) 



1/2 



(7) 



(8) 



(9) 



with 



^=Jl+[2 V /r V'(r )Y 



On the other hand, vq is chosen to minimize Eq, such that, 



dE n , d 2 E 



dro 



= and 



drl 



> 



and then one can get 

1 + 2 £ + M (2n + l)u = r 2 V'(r ) + 2 AA) 1 '\ 



(10) 



(11) 



(12) 



which is an explicit equation in r$. It is convenient to summarize the above procedure in 
the following steps: (a) Calculate r Q from Eq.(12) and substitute it in Eq.(9) to find Q 
(b) Substitute Q in Eq.(3) to obtain E Q . (c) Finally, one can obtain E Q and then calculate 
E n i from Eq.(6). However, one is not always able to calculate r Q in terms of the potential 
coupling constants since the analytical expressions become algebraically complicated, al- 
though straightforward. Therefore, one has to appeal to numerical computations to find 
r D , and hence E Q . 



3 Application, Results and Discussion 



The dilaton is a scalar field predicted by superstring theory [12]. The mechanism and the 
form of the dilaton potential are unknown, although it is believed that they could be related 
to nonperturbative sector of the theory. Recently it was observed in [5] (see also [6,7]) that 
a string inspired coupling of a dilaton cj) to 4d SU{N C ) gauge fields = T a A a , with T a 
the (AT? — 1) SU(N C ) generators, yields a phenomenologically interesting potential V(r) 
for the quark-antiquark interactions. Dick interquark potential was obtained as follow: 
First start from the following effective field theory with the dilaton-gluon coupling G((f>) 
and the dilaton potential W ((/)): 

= W - \(d^? + W(d>) + J$A* (13) 
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then construct G(<f) under the requirement that the Coulomb problem still possesses 
analytical solutions. The coupling G{(j>) and the potential W((/>) that emerged are: 

f 2 1 
GU) = const. + J -r, W(<f>) = -m 2 <p 2 (14) 
qr 2 

where / is a scale parameter characterizing the strength of the scalar-gluon coupling and 
m represents the dilaton mass. 

Next, consider the equations of motion of the fields and 4> and solve them for 
static point like color source described by the current density J£ = p a rj^°- After some 
straightforward algebra, Dick shows that the interquark potential Vo(r) is given by (up 
to a color factor), 

V D (r) = +5/ yZ^Z ln[exp(2mr) _ 1] (15) 

Eq.(15) is remarkable since at large values of r it leads to a linear confining potential 
Vo{r) ~ 2gfm^ ^_^ r. This derivation provides a challenge to monopole condensations 
as a new quark confinement scenario. Therefore, it is well justified to dedicate more efforts 
to the investigation of this effective coupling function G(4>) and to the phenomenological 
application of Dick potential Vo(r). The method of obtaining results from the theory 
requires us to choose several numerical inputs, examination of the previous section shows 
that there are five parameters to be calculated, namely, m c , rrib, m, f and a s . A few 
comments about the choice of the parameters are in order. The numerical values of the 
charmed-quark mass m c = 1.89 GeV and the bottom-quark mass nib = 5.19 GeV are 
chosen so that, the QCD coupling constant a s (fi) at any renormalization scale can be 
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calculated from the world average experimental value a s (m z ) = 0.117 via [13] 



a s (m z ) 



(16) 



1 - (11 - \n f )[a s (rn z )/2ir\ln(rn z /iiy 



then one has 



a s (m c ) = 0.31, 



a s (m b ) = 0.2. 



(17) 



For the be quarkonia, we obtain a s (4^ c ) = 0.22, where fi bc is the reduced mass [14]: 



On the other hand, the potential parameters m and / are considered free in our analysis 
and are obtained by fitting the spin-averaged cc, and bb mesons. An excellent fit with the 
available experimental data can be seen to emerge when the following values are assigned 



The results of our calculation for the spin-averaged energy levels of interest are given in 
Tables (1,2). In all cases, where comparison with experiment is possible, agreement is very 
good. We also present, in Table 3, the results for the be quarkonia. Our estimate for the B c 
mass, the lowest pseudoscalar S-state of the spectra, is compatible with the experimental 
value reported in [15]. As to the higher states masses, they compare favorably with other 
predictions based on QCD sum-rules [16] or other potential models [17]. In conclusion, 
to the best of our knowledge, this is the first time that Dick interquark potential of 
Eq.(15), is tested and used successfully to fit the spin-averaged cc, bb, and be systems. 
This potential will certainly open a new window in the potential model applications, 



m b m c 



(18) 



Hbc = 



77i5 + m c ' 



m = 56.9 MeV, 




(19) 
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essentially those with a Coulomb plus a confining potential. On the other hand, since a 
unique theory/scenario for the dilaton mass is still lacking [18], our estimate for the mass 
of the dilaton m = 56.9 MeV, resulting from a fit to existing experimental data of heavy 
quarkonium, lies in the range given in [19]. As suggested by the authors of Ref.[20], the 
possibility to identify the dilaton particle to a fundamental scalar, invisible to present day 
experiment, should not be ruled out. 
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Table 1: Calculated spin-averaged data (in units of GeV) M n £ of charmonium cc energy 
levels in the semi-relativistic wave equation using Dick potential 



State, n£ 


M n£ , SLET 


M n£ , Exp. 


State, ni 


Mru, SLET 


M ni , Exp. 


IS 


3.073 


3.068 


IP 


3.548 


3.525 


2S 


3.662 


3.663 


2P 


3.871 




3S 


4.027 


4.028 


ID 


3.787 


3.788 



Table 2: Calculated spin-averaged data (in units of GeV) M n i of bottomonium bb energy 
levels in the semi-relativistic wave equation using Dick potential 



State,77i 


M n£ , SLET 


M n £, Exp. 


State,rii 


M ne , SLET 


M n i, Exp. 


IS 


9.450 


9.446 


IP 


9.903 


9.900 


2S 


10.014 


10.013 


2P 


10.206 


10.260 


3S 


10.292 


10.348 


ID 


10.129 





Table 3: Calculated spin-averaged data (in units of GeV) M n £ of be energy levels in the 
semi-relativistic wave equation using Dick potential 



State,n£ 


M n£ , SLET 


M n£ , Exp. 


State,ra£ 


M ne , SLET 


M n£ , Exp. 


IS 


6.322 


6.40 ±0.39±0.19 


IP 


6.767 




2S 


6.876 




2P 


7.072 




3S 


7.161 




ID 


6.994 
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